First-Principles Modeling of Pt/LaAlOs/SrTiOs Capacitors 
Under an External Bias Potential 
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We study the electrical properties of Pt /LaAlOs / SrTiOa capacitors under the action of an external 
bias potential, using first-principles simulations performed at constrained electric displacement field. 
A complete set of band diagrams, together with the relevant electrical characteristics (capacitance 
and built-in fields) , are determined as a function of LaAlOa thickness and the applied potential. We 
find that the internal field in LaAlOa monotonically decreases with increasing thickness; hence, the 
occurrence of spontaneous Zener tunneling is ruled out in this system. We discuss the implications 
of our results in the light of recent experimental observations on biased LaAlOa/SrTiOs junctions 
involving metallic top electrodes. 

PACS numbers: 73.20.-r, 71.15.-m, 71.70.-d 



In recent years oxide-oxide hetero junctions have gen- 
erated widespread interest, because of their potential 
for applications in microelectronics and intriguing funda- 
mental physics. When thin films of polar LaAlOa (LAO), 
composed of formally charged (LaO)+ and (A102)~ lay- 
ers, are stacked on top of the non-polar (001) surface of 
SrTiOa (STO), a two-dimensional metallic electron gas 
appears at the interface under suitable experimental con- 
ditions [1^. In particular, Ti02-terminated LAO/STO in- 
terfaces are found to undergo an insulating-to-conducting 
transition (named "electronic reconstruction" ) when the 
thickness of the LAO film exceeds 3 — 4 layers Al- 
ternatively, 2D metallicity can be switched on and off by 
applying an external bias potential between the conduct- 
ing LAO/STO interface and a metallic electrode placed 
on top of the free LAO surface [3-8]. In the latter setup 
(i.e. an asymmetric capacitor where the bottom electrode 
is the 2D electron gas, see Fig. [T]) novel striking phenom- 
ena have been observed recently, which hold promise for 
the realization of novel field-effect devices with high op- 
eration speed and low power consumption. 

First, Singh-Bhalla et al. Q reported a non-trivial de- 
pendence of the tunneling current on the LAO thickness, 
^LAOi with an abrupt increase at c^lao — 20 unit cells. 
This was ascribed to the presence of a built-in electric 
field, which would cause Zener breakdown in the LAO 
film after exceeding a certain c^lao critical value. The 
electrical response of the tunnel junction also displayed a 
clear hysteretic behavior as a function of applied voltage, 
suggesting interesting opportunities for memory applica- 
tions. Second, the measured capacitance was found to 
undergo a remarkable increase at low carrier densities Q . 
This enhancement was ascribed to the peculiar electronic 
properties of the 2D electron gas at the LAO/STO inter- 
face, and in particular to a "negative compressibility" 
regime. While such a behavior was already known in the 
context of semiconductor heterojunctions, its unprece- 
dented magnitude (40%) in the LAO/STO system chal- 
lenges the current theoretical understanding of this effect. 




£(STO) 



£ d 

LAO LflO 



FIG. 1: (a) Schematic illustration of the Pt/LAO/STO ca- 
pacitor; large blue spheres represent Sr atoms, small blue 
Ti, red O, large purple La, pink Al and small purple Pt. 
(b) Pt/LAO/STO band diagram; (^i^'^' are the n-type SBH 
at the LAO/STO and Pt/LAO junctions, £lao is the electric 
field in LAO, dLAO is the thickness of the film, Ey / Ec are 
the valence / conduction band edges, Ef are the Fermi levels, 
and Vcxt is the external bias potential. 



Rationalizing these phenomena in terms of the micro- 
scopic properties of the LAO/STO and LAO/electrode 
junctions is very desirable in sight of future progress. In 
this context, first-principles electronic structure methods 
appear ideally suited to describing, with unbiased accu- 
racy, the subtle interplay between carrier confinement, 
polar distortions and electrical perturbations applied to 
the sample. Indeed the LAO/STO system has been ad- 
dressed by a large number of ab-initio studies in the past 
few years [9 16]. However, studying the phenomena de- 
scribed in Refs. [7| and Q entails some additional tech- 
nical challenges, because of the necessity of introducing 
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FIG. 2: (a) Layer-resolved La 5s semicore energies, -Bj^'^', 
obtained at D = -0.5 e/5 (left) and -0.3 ejS (right). Trian- 
gles represent the extrapolated values at the interfaces, 4>tA^ ■ 
(b) Sketch of the Pt/LAO/STO band diagram and details of 
the extrapolation scheme (see text). 



an external voltage in a complex capacitor system that 
is overall metallic. 

In this Letter, we propose an elegant solution to this 



problem by using recently developed methods [17|, |18| 
to perform ah initio simulations at constant electric dis- 
placement D. In particular, we show that: (i) the in- 
ternal electric field (^lao) in a Pt/LAO/STO capaci- 
tor is not an intrinsic and fixed quantity, but a func- 
tion of both applied bias and LAO thickness; (ii) at zero 
bias, ^LAO monotonically decreases with film thickness as 
~ l/dLAO, so that no Zener breakdown ever occurs; (iii) 
the LAO/STO interface is associated with a capacitance, 
predominantly due to band bending effects, which is con- 
stant within a wide range of carrier densities. Based on 
these findings we deduce a complete set of ah initio band 
diagrams and a simple analytical expression for f lao ^ 
a function of the relevant parameters. Finally, we dis- 
cuss the implications of our results for the interpretation 
of experimental observations in [§] and ■ 

To start with, we are interested in determining the 
electronic and structural properties of an arbitrarily thick 
Pt/LAO/STO capacitor subjected to an external bias po- 
tential Voxt • Instead of working directly at fixed T4xt , we 
use the electric displacement field D as independent elec- 
trical variable; this has clear advantages from the point 
of view of modeling [l^ , as it allows one to break down a 



layered system into smaller constituents and treat them 
separately ("locality principle"). (Note that working at 
fixed V or at fixed D yields the exact same informa- 
tion, except that we express it in terms of two different 
macroscopic electrical variables; these are related by a 
Legendre transformation [l7|.) In our case, we shall con- 
sider the LAO/STO interface [(1) henceforth], bulk LAO 
and the Pt/LAO interface [(2) henceforth]. Following the 
arguments of Refs. [13] and [23], we can decompose the 
total potential drop Vtot across the capacitor into three 
terms (see Fig. [1]), 



V,ot{D,N) - <i>^^^{D) - ^\tHD) + NVi^Ao{D) . (1) 

Here T^lao = ^ i^^lao clao (clao is the out-of-plane lat- 
tice parameter) is the potential drop across one unit cell 

(1 2) 

of bulk LAO, TV the number of unit cells, and (f>n ' 
are the (r'-dependent) n-type Schottky barrier heights 
(SBH) at the metal/insulator junctions 17[. Note that 



to obtain the ground state of the system at a given ap- 
plied potential Vcxt one simply needs to invert Eq. ([1]) 
and solve for Vtot{D,N) = Vcxt- This way, the daunt- 
ing problem of simulating the full Pt/LAO/STO capac- 
itor under an external bias reduces to the more familiar 
task of calculating SBH at metal/insulator interfaces as 
a function of D. As we shall explain in the following, it 
is relatively straightforward to do this with a standard 
first-principles code, without the need for a specialized 
finite-field (or even Berry-phase) implementation. 

Our calculations are performed within the local den- 
sity approximation of density functional theory and the 
PAW method 2l|, as implemented in the in- house code 
Lautrec [i^]. We compute the quantities on the right 
hand side of Eq. ([1]) by means of three separate calcu- 
lations: two interface models within a X/LAO/vacuum 
slab geometry, where X is the metallic electrode (either 
Pt or STO), and a periodic LAO bulk model. For the 
X/LAO interfaces we use stacks of 10/7 (X=Pt) and 
12/5 (X=STO) layers respectively. The in-plane peri- 
odicity is set to 1 X 1 perovskite cells, with the lattice 
parameter fixed to the theoretical equilibrium value of 
bulk STO (asTO = 3.85 A). To constrain the electric dis- 
placement to a given value we introduce a layer of bound 
charges Q at the free LAO surface via the virtual crystal 
approximation 23]. By applying a dipole correction in 
vacuum, we enforce D — outside the free surfaces; then, 
provided that the surfaces remain locally insulating, we 
have D = Q/S We explore values of D within 

the range —0.5 e/S < -Dlao < —0.3 e/S (Dsto is set to 
zero). 

To extract the dependence of the SBH on D, we use 
computational techniques similar to those in Refs. [i^ 
and [24]. In particular, for the estimation of the band off- 

(1 2) 

set (j)n ' (D) we need two independent quantities, the so- 
called lineup term (an interface property) and the band- 
structure term (a bulk property) [24| . The lineup term 
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FIG. 3: (a) Calculated n-SBH as a function of D. The curves 
represent the linear fits (t)\t\D) = 0.35 -I- 1.62 (D -I- e/2S) 
and (f>'ii\D) = 1.18 + 0.05 (i3 + e/2S) (in units of V). 
(b) Calculated potential drop per unit cell of bulk LAO; 
VLao(^) = 1-27 {D + e/2S) + 2.24 (D + e/2Sf (solid line) 
and Vi,Ko[D) = 1.36 (D + e/2S) (dashed line), (c) ^lao as 
a function of the number of LAO layers and bias potential 
—0.5 < l/ext < 0.5 V. The i-axis is in logarithmic scale. Sym- 
bols and lines represent results obtained with Eqs. ^ and ((6]), 
respectively. Inset: Static dielectric constant of bulk LAO as 
a function of DS/e (solid curve) and the corresponding aver- 
aged value £lao (dashed line). 



relies on the choice of a reference energy in the insulating 
LAO layer; here we use the La 5s semicore energy, that 
we extract from the layer-resolved density of states of the 

(1 2) 

slab models. We define as E^^' {D,j) the energy loca- 
tion of the La 5s peak in the j-th LaO layer of the (1) 
or (2) system, referred to the respective Fermi level. In 
the case of zero field, £^5^,' {D,j) converges to a constant 
value for layers j lying sufficiently far away from the in- 
terface or the surface (typically two-three unit cells), and 
the definition of the lineup term is straightforward. Con- 
versely, when the electric field in LAO is nonzero (see 
Fig. [5]), E^s has a linear dependence on the layer index. 
In particular, deep enough in the film, the variation of 
i?5s between two consecutive cells corresponds precisely 
to the bulk internal field at the same value of D, 



E. 



(1,2)/ 
5s 



iD,j + l)-Et'\D,j)^V^AoiD). (2) 



Therefore, once Vlao(-D) is known (from the bulk calcu- 
lations), we can define an extrapolated reference energy 
(see Fig. [2]), 



El'/\D,j)^k{D)+jVi^Ao{D), 



(3) 



where the constant k{D) is chosen in such a way that 
E,-ci{D,j) — E5s{D,j) far from the interface. Then, we 
define the lineup term as 



where ' ' indicates the interface plane location; here we 
assume j| = 3/4 and = 1/4 (see Fig.©. Finally, to 
obtain the n-type SBH </>„ ' (D) we just need to add to 
'^rof (^) the D-dependent band-structure term [we cal- 
culate it in the bulk, as the difference between the La(5s) 
level and the conduction band minimum], -Ecbm(^), 



<t^i''^HD)^^ilfHD)+EcBM{D). 



(5) 



f 1 2) 

Our results for 0„ {D) and Vlao(-D) are shown in 
Fig.|3lja,b). We find that (f)n'^\D) behave linearly within 
the whole studied range. Consequently, the total po- 
tential drop due to both interfaces can be expressed as 
Vi{D) = + Cfi iD + e/2S), where is the total 
interface potential drop at I? = —e/2S and the coeffi- 
cient Cj"^ physically corresponds to the overall inverse 
interface capacitance density 17, [2^. Considering the 



potential drop in bulk LAO also to be linear in D [see 
Fig. [Sib)], i.e. Flao(-D) ~ {D + e/2S), we obtain 
a simple yet very accurate analytical expression for the 
electric displacement field in LAO, 



D{V,.t,N) + — 



(6) 



(1,2), 



,-(l,2)> 



(4) 



where C^^ = 1.36 mVF, = -0.83 V, and Cf ^ = 
1.57 m^/F {S = flsTo)- The inverse bulk LAO capaci- 
tance density, C^^ , is directly related to the static dielec- 
tric constant of LAO, elao, via C^^ = clao / (eoCLAo) ; 
in our computational model we have elao = 31. For a 
given value of _D, we compute the corresponding electric 
field using ^lao = (-D + e/2S) /elao- [In the "exact" 
treatment we replace the average clao with the calcu- 
lated eLAo(^), see the inset of Fig. WL^)]- Our results 
for £LAo(K:xt, are shown in Fig. [3][c). It is worth 
noting that the outcomes of Eq. ^ [exact, considering 
non-linearities in Vlao(-C)] and Eq. ([6]) (approximate) 
are in excellent agreement. In both cases we find that 
^hAoiVc^it, N) monotonically decreases with the inverse 
of LAO thickness, independently of the applied bias po- 
tential. 

These results have important implications concerning 
the interpretation of the experiments reported in Ref. ■ 
First, no intrinsic built-in LAO electric field exists in the 
short-circuited Pt/LAO/STO capacitor system, contrary 
to the assumptions of Singh-Bhalla et al. [7| . Second, and 
most importantly, the uniform decrease of £'lao {Vcy^t , N) 
with LAO thickness in practice rules out the hypothesis 
of Zener tunneling that was used in Ref. [3] to explain 
the abrupt increase of the tunneling current at A'^ 20. 
Indeed, as fLAo(K:xt, decreases with increasing iV, 
Zener tunneling becomes increasingly less likely in thicker 
films; note that in our model the LAO valence band never 
goes above the STO Fermi level unless a strong external 
potential Kxt ^ 1.8 — 1.9 eV is applied. An alternative in- 
terpretation of the aforementioned tunneling experiments 
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FIG. 4: (a) Mesoscopic band diagram of the LAO/STO 
junction. The electrostatic potential in STO, Vh, (thick 
solid curve), the conduction and valence band edges (shaded 
dashed curves) and the Fermi level (thick dashed line) are 
shown, (b) Ti t2g density of states of our TB model {Vn is 
constant and set to zero), (c) Vh as calculated with the TB 
model (Fermi levels are shifted to zero), (d) Electrostatic po- 
tential at the interfacial Ti site as a function of D. Empty 
circles (filled squares) are calculated for a 12 (24)-layer STO 
slab and the solid line represents a linear fit to the data. 



rests on the results of Ref. [26| . Precisely at a thickness 
of 20 unit cell, Cancellieri et al. reported the onset 
of in-plane strain relaxation in the LAO overlayer. Such 
relaxation processes are known to induce defec ts ( e.g. 



cracks, or misfit dislocations) in the LAO film [26|, l27|. It 
is not unreasonable to think that these defects might con- 
stitute preferential paths for electron or hole conduction; 
also, these defect-mediated conduction processes might 
explain the hysteretic behavior of the electrical diagrams 
of Ref. |3] . (The authors invoked the presence of switch- 
able dipoles at the LAO/STO interface, which appears 
puzzling as neither LAO nor STO are ferroelectric.) 

Finally, we shall comment on the microscopic mech- 
anisms that lead to the calculated interface properties. 
The total interface capacitance density can be expressed 
as the sum of two contributions, 



^-1 _ d 



(2)^ 



1 



1 



C(i) C(2) 



(7) 



(2) 

each due to a different interface. The variation of <f>n {D) 
with D is very small (see Fig. [3]); this means that the Pt 



electrode behaves ideally ("perfect" screening [2J|), i.e. 
1/C(2) ^ 0, and therefore we restrict the following anal- 
ysis to the LAO/STO junction. Specifically, we shall 
consider the effects of band-bending on C^^) , i-e. the D- 
dependent interfacial potential drop that occurs in STO 
due to the presence of the confined charge carriers [l^l ■ 
To estimate the band-bending contribution, (^sto(-P), we 
have employed the tight-binding (TB) model |29[ intro- 
duced in Ref. [18]. We define 0sTo(-D) as the value of 



the electrostatic potential in STO at the interfacial Ti 
site, referred to the Fermi level of the confined electron 
gas, at a given D [see Fig. |4l[a,c)]. Our results show that 
4'STo{D) behaves linearly within a wide range of val- 
ues [Fig. mjd)]. The band-bending contribution to Cj^^ 
corresponds to the slope of the linear fit, which amounts 
to 1.03 m^/F. This contribution constitutes a significant 
fraction (approximately two thirds) of C'^^. The remain- 
der of Cj^J is likely to be caused by other effects that 
typically occur at oxide-oxide interfaces, e.g. intrinsic 
dielectric "dead layers" S^, Slj. As these appear to be 
of secondary importance here, we have not pursued this 
analysis further. 

The above results are timely for the interpretation 
of recent capacitance measurements of thin-film elec- 
trode/LAO/STO heterostructures Q. In the regime of 
large carrier densities, Li et al. measured a dielectric 
constant of the LAO film of eLAO ^ 18, which is signif- 
icantly smaller than the typical value reported for LAO 
single crystals (eLAO = 25 — 30). Our results indicate 
that this discrepancy may be related to interfacial ef- 
fects, and in particular to the relatively small value of 
C*-^-', rather than to the quality of the LAO film. It is 
worth noting that Li et al. [8j] used YBCO electrodes, in- 
stead of the Pt electrodes considered here. This implies 
that a possible contribution from the upper electrode in- 
terface should be in principle taken into account (recall 
that 1/C*^^-' essentially vanishes in the present Pt/LAO 
case) when comparing our data to those of Ref. |8|. In 
this context, note that Eq. ([6]) is completely general and 
applies to any electrode/LAO/STO configuration; one 
simply needs to replace two parameters (that are specific 
to the electrode/LAO interface), the Schottky barrier at 
zero field (jy^n^ and the interface capacitance C^^-' . In other 
words, one does not need to repeat the calculations of the 
LAO/STO interface, since we have essentially decoupled 
the bottom from the top electrode and their properties 
can be computed separately. 



In summary, we have studied the electrical properties 
of metal/LAO/STO capacitors fully from first-principles, 
determining a complete physical picture of the band off- 
sets and internal fields. Our results and methodologies 
open new avenues in the first-principles study of func- 
tional oxide heterostructures, and provide useful guide- 
lines for the interpretation of the available experimental 
data on this system. 
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